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Wingbeat  Forcing  Functions:  Part  II  -  Controller 
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A  control  strategy  is  proposed  for  a  minimally-actuated  flapping-wing 
micro  air- vehicle  (FWMAV).  The  proposed  vehicle  is  similar  to  the  Har¬ 
vard  RoboFly  that  accomplished  the  first  takeoff  of  an  insect  scale  flapping 
wing  aircraft,  except  that  it  is  equipped  with  independently  actuated  wings. 
Using  the  derivation  of  the  aerodynamic  forces  and  moments  from  Part  I, 
a  control  allocation  strategy  and  a  feedback  control  law  are  designed  that 
enable  the  vehicle  to  achieve  untethered,  stabilized  flight  about  a  hover 
condition.  Six  degree-of-freedom  maneuvers  near  hover  are  demonstrated 
as  well.  The  control  laws  are  designed  to  make  use  of  two  actuators  that 
control  the  angular  position  of  the  wing  in  the  stroke  plane.  The  Split- 
Cycle  Constant-Period  Frequency  Modulation  with  Wing  Bias  technique, 
introduced  in  Part  I,  is  used  to  allow  each  wing  to  generate  non-zero  cycle- 
averaged  aerodynamic  forces  and  moments.  This  technique  modifies  the 
frequencies  of  the  up  and  down  strokes  to  yield  non-zero  cycle- averaged 
drag  due  to  the  flapping  motion  of  a  wing.  Additionally,  the  midpoint  of 
the  wingbeat  profile  can  be  modified  by  use  of  a  wing  bias.  The  bias  is 
introduced  to  primarily  provide  pitching  moment  control.  In  this  work, 
the  sensitivities  of  cycle-averaged  forces  and  moments  with  respect  to  the 
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control  parameters  are  computed.  Single  axis  controllers  are  designed  and 
the  complete  system  is  simulated. 

I.  Introduction 

In  Part  I1  of  this  series,  expressions  for  the  instantaneous  and  cycle-averaged  aerodynamic 
forces  and  moments  of  a  flapping  wing  biomimetic  vehicle  were  developed.  This  vehicle  is 
similar  to  the  Harvard  Robofly  described  by  Wood,2  with  the  primary  exception  being  the 
proposed  vehicle  makes  use  of  two  actuators,  which  are  used  to  independently  vary  the 
position  of  the  wing  spars  in  the  stroke  plane.  In  previous  work,3,4  a  third  actuator  was  used 
to  control  the  position  of  a  bobweight  such  that  the  center-of- gravity  could  be  manipulated 
to  control  the  pitching  moment.  In  this  work,  the  third  actuator  has  been  eliminated  by  the 
inclusion  of  a  wing  bias  term.  The  goal  of  this  paper  is  to  evaluate  the  cycle-averaged  control 
derivatives  with  respect  to  all  control  variables  to  determine  if  sufficient  control  authority  can 
be  achieved  to  regulate  the  vehicle’s  6  degree-of-freedom  (DOF)  body  position  and  attitude 
using  wingbeat  frequencies,  curw  and  ojlw,  the  split-cycle  parameters,  8rw  and  8lw,  and 
the  wing  bias  terms,  t]rw  and  t]lwi  as  control  inputs. 

In  the  present  paper,  the  analysis  that  was  started  in  Part  I1  is  continued  by  analyzing  the 
effects  of  wing  bias  and  other  control  parameters  upon  the  body-axis  forces  and  moments. 
Expressions  for  the  cycle-averaged  control  derivatives  associated  with  the  variations  in  the 
fundamental  wingbeat  frequencies,  split-cycle  parameters,  and  wing  biases  are  derived.  The 
expressions  are  linearized  about  a  hover  flight  condition  and  a  control  allocation  strategy 
is  proposed.  A  controller  is  developed  such  that  the  vehicle  can  perform  waypoint  tracking 
and  simulation  results  are  presented. 

II.  Control  Derivatives 

The  parameters  used  to  control  the  aerodynamic  forces  and  moments  are  the  fundamental 
wingbeat  frequencies,  ujrw  and  colw,  the  split-cycle  parameters,  8rw  and  8lw,  and  the  wing 
bias  terms,  r]RW  and  r/LW.  For  controllability  analysis  and  control  synthesis,  the  sensitivity  of 
each  cycle-averaged  force  and  moment  to  each  control  input  parameter  must  be  determined. 
The  forces  and  moments  are,  in  general,  functions  of  the  control  input  parameters  as  well  as 
the  past  values  of  r]RWl  r]LW  because  when  ARW  =  ALW  =  1,  A At  =  rji  —  rj^ i  =  77  —  77-1.  Note 
that  the  subscript  i  notation  has  been  removed  and  it  should  be  understood  that  i]  =  77*  (the 
current  value  of  wing  bias)  and  77-1  =  rji-i  (the  one  cycle  past  value  of  wing  bias).  Hence, 
the  cycle-averaged  forces  and  moments  are  sensitive  to  changes  in  77_1?  although  77^  is  not 
considered  a  control  input.  Letting  G  be  a  generalized  cycle- averaged  force  or  moment,  the 
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total  increment  of  G  is 
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A  feedback  controller  will  generate  a  set  of  desired  cycle-averaged  forces  and  moments,  hence, 
the  left-hand  side  of  Equation  1  will  be  replaced  A Gdes-  Of  interest  is  the  hover  condition 
and  Equation  1  can  be  written  as 
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The  increments  in  some  parameters  are  replaced  by  the  parameters  themselves  since,  at 
hover,  5RW  =  5lw  =  Vrw  =  Vlw  =  Vrw- i  =  Vlw- i  =  0.  The  partial  derivatives  in 
Equation  2  are  now  evaluated.  The  expressions  for  the  cycle-averaged  forces  and  moments 
are  given  in  Part  I1  of  this  series. 


A.  X-Body  Axis  Force  Control  Derivatives 
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For  the  x-body  axis  force  sensitivity  with  respect  to  Vrw,  recall  that  A Arw  is  a  function  of 
Vrw-  Therefore, 
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Left  wing  expressions  are  of  the  same  form  as  the  right  wing. 
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B.  Y-Body  Axis  Force  Control  Derivatives 
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(14) 


where  Ji(-)  is  a  Bessel  function  of  the  first  kind  and  H i(-)  is  a  Struve  function.5 
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C.  Z-Body  Axis  Force  Control  Derivatives 
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+  {Arw  +  A  Arw)  - — — — — tt - cos  rjRW  —  Ji{Arw  +  A  ARW)  sin  r]RW 


dF. 


-B 

zRW 


9Hi{ArW  +  A  Arw) 
drjRw 

kn^'RW  {urw  +  cfrw ) 


dljRW-l 


dr]Rw 

si  v-Vrw  ~  Hi{Arw  +  A  Arw)  cos  rjRW 


Ji{Arw  +  AArW)  COS  T)rw 


(21) 


TT  /  A  »  a  \  ■  l  9AArW 

—  Hi  {Arw  +  A  Arw)  sin  r/Rw  \  - 

J  9t)rw- i 

!  A  A  A  I  f  9Ji{Arw  +  AArw) 

+  {Arw  +  A  Arw)  < - — - - cos  rjRW 


9Hi{Arw  +  A  Arw) 

9t]rw-  1 


sm  tjrw 


(22) 
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d  F. 


-B 

zLW 


knAjwh  {Alw)^lw  cos  Vlw  kp 


d5 


LW 


U 


LW 


4  \(ujLw  ~  ^lwY 

*  [Ji{Alw)  cos tilw  +  Hi(Alw)  sin ijLw] 

+  {Alw  +  A Alw)  J\(Alw  +  A Alw)  cos t]lw 
—  Hi  ( Alw  +  A Alw)  sin  rjLW 


(23) 


OF- 


-B 

zLW 


duj 


LW 


kpALWJi(ALW)  (2a jLW  —  §LW )  cos rjLw  kn  ( 2c o\w  —  75lw^lw  +  4a Jlw^lw 

2  4  \  (lolw  ~  25 lw )2 


dF- 


-B 

zLW 


*  i  Alw  [Ji(Alw)  cos  ijLw  +  H\{Alw)  sin t]lw\ 
+  {Alw  +  A  Alw)  J\{Alw  +  A  Alw)  cos  rjLw 

—  Fti  {Alw  +  A  Alw)  sin  rjLw 
—kpALw  Ji{ALw)ulw  {^lw  ~  &lw)  sinrjLw 


(24) 


drjLW 


kpALW^LW  (k ’LW  +  aLw) 


—  J\{Alw )  sin  r/Lw  +  Hi{ALw )  cos t/lw] 


kp^LW  {^LW  +  &Lw) 


H\{Alw  +  A  Alw)  sin  tjlw 


J\{Alw  +  A  Alw)  cos  t/lw 
]  dA  ALw 


j  dr]Lw 


+  {Alw  +  AAlw)  ^  ^ — ~rr - cos  t/lw  ~  J\{Alw  +  A ALw)  sin  rjLw 


dF. 


zLW 


drjLW- i 


dljLW 

9Hi{ALw  +  AALw)  .  u  ,  A  ,  A  A  ^ 

- - - sin  rjLw  ~  Hi(ALw  +  AALw)  cos  t/lw 

dr/Lw 

kpUJLW  {^LW  +  &lw)  I  j  J  ,  A  A  a  I 

Ji{Alw  +  A  Alw)  cos  ijLw 
)  8A ALw 


(25) 


H\{ALw  +  A ALw)  sin rjLw 


J  drjLW -l 


,  ,  A  ,  A  \  i  )  dJ\{ALW  +  AAlw) 

+  {Alw  +  AAlw)  ^ - b - cos  rjLw 

orjLW-i 


dHi{ALw  +  A  ALw) 
drjLW -i 


sin  rjLw 


(26) 
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D.  Rolling  Moment  Control  Derivatives 

The  rolling  moment  control  derivatives  are 


9MxRW  kD^-RW^RW  D .  ,T 

ft  r  —  ~a  Kl\ll RW 

OOrw  4 

— ^  [Arw RM'2rw  +  (ArW  +  A Arw)  RM3rw\  - HA. — 

°  \Urw  —  2orw 


where 


(27) 


RMirw  =  yZP Arw  +  J\  {Arw)  { w  cos  r)RW  +  2 A zR  sin  rjRW  }  (28) 

RA12rw  =  vZ/Arw  +  w  {  Ji(ArW)  costirw  +  Hi(Arw)  sin rjRW} 

+2A zR  {J\{Arw)  smrjRw  —  Hi{Arw )  cos  i]rw}  (29) 

RM3  rw  =  yZP  ( ArW  +  AArw)  +w  |  J\  ( Arw  +  A  Arw)  cos  r]RW 

—Hi  {Arw  +  A  Arw)  sin  tjrw  1 

+2A  zR  {Ji  {Arw  +  AArw)  sin  t)rw  +  Hi  {Arw  +  A  Arw)  cos  t/rw]  (30) 
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dr/Rw 


kDARW( ^nw  -  inw)  RMikw  _  kn  [ArwRm^  +  (Amr  +  AArw)  rm3rw) 

:|.  —  76rw^rw  +  4c orw&rw) 

{urW  —  2  5rw) 

koARwURW  (urw  -  Srw)  r  J  ,  A  .  ,oaB  'll 

- - - [J\{Arw)  { -w  sm  r/Rw  +  2 A zR  cos  yRw)\ 

^dArW^RW  (^RW  +  (Jrw)  r  T  t  A  \  ■  \  TT  (  A  \  1 

- - -  W  {  —  Ji{ArW)  Sill  yRW  +  Hi(ArW)  COS  1]RWJ 


2A ZR  {J\{Arw)  cos  i]rw  +  Hi(Arw)  sin  t/rw } 
kp^Rw  {urW  +  cr rw)  dAARW  \wp 


dijRW 


y cp  ( Arw  +  aarw) 


+  w  {Ji  ( ArW  +  A Arw)  cos  r]RW  —  Hi  ( Arw  +  A Arw)  sin  t]rw } 

kp^RW  (URW  +  &Rw)  /  ,  A  A  \\  WP^^Arw 


(Arw  +  AArw)  ycp 


dyRW 


f  8Ji(Ariy  +  A  Arw) 

+W  | - — - COS  7j RW  —  J\  ( ArW  +  AArW)  Sin  1]RW 


3Hi(ArW  +  AArw) 


'  inw 

OT]rw 

_kDujRw  (‘ aJrw  +  aRw)  8AArw 

8  dr]RW 

+  Hi  ( Arw  +  A Arw)  cos  t]rw  ] 


sin  Vrw  ~  Hi  ( ArW  +  AArW)  cos  t/rw 


2Azft  <  Ji  (Arw  +  AArW)  si nrlR\v 


kp^RW  (URW  +  &Rw)  /  ,  ,  A  A  \  [oA  B  f  9Ji(Arw  +  AArw) 


(Arw  +  AArw)  2A zr 


dyRW 


sm  yRW 


,  .  .  .  ,  dHi(ARW  +  AArw) 

+J\  (Arw  +  AArw)  cos  yRw  H - tt - cos  t/rw 

oriRw 

—  Hi  (Arw  +  AArw)  sin  t/rw  1 
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-B 


9MxRW  kjjujRw  (ujrw  +  & rw)  9AARW 


drjRw- 1 


driRw- 1 


2/™P  (^4-Rw  +  AArW) 


yep 


+iu  {  J\  ( Arw  +  AARW)  cos  tjrw  —  Hi  (. ARW  +  AARW)  sin  t]rw} 


^D^RW  (Wrw  +  Orw) 


{Arw  +  AArw) 


wp 
t/  cp 


OAA 


RW 


drjRw- 1 


w 


9Ji(Arw  +  AArw ) 


driRw -1 

—  knujRw  {urw  +  ctrvk)  9AArw 
8  drjRw- 1 

+  #1  (^4rw  +  AARM/)  cos  rjRW 

^d^rw  {urw  +  ctrw) 


cos  rjRW  - 


dH\  (  Arw  +  AARiy) 


9vrw-  1 


sm  rjRW 


2A 2^  Ji  (Arw  +  AArW)  smrjRW 


{Arw  +  AArw) 


2Az 


R 


9Ji(Arw  +  AA/jvn) 
9tjrw-  1 


9H\  (  Arw  +  AA/jvk) 

H - o - cos  r]RW 

9t)rw-  1 

For  the  left  wing,  the  rolling  moments  control  derivatives  are 

_ q 

9MX  LW  kc>ALWujLW 


■  sm  r]RW 
(33) 


LW 


A: 


+  wr  [AlwRM2lw  +  +  AALw)  RM3lv 


CO 


(34) 


nr 


(wnr  —  ^lw) 


where 


RMilw  =  yWpALW  +  J\{Alw)  [w cos riLW  +  2Azf  sin  r/LH/ }  (35) 

RHhLW  =  y^vF Alw  +  w  {  Ji(Alw)  cos t]lw  +  Hi(Alw)  sinr^vv} 

+2Azf  {J\(Alw)  sin  i]LW  —  Hi(Alw)  cos  t/lw}  (36) 

RM3lw  =  yWp  ( Alw  +  A Alw)  +  w  |  J\  (Alw  +  AAlw)  cos  r]LW 

— H\  ( Alw  +  A Alw)  sin  r]LW  | 

+2A zR  {  J\  ( Alw  +  A Alw)  sin  yLW  +  Hi  ( ALW  +  AALW)  cos  ijlw}  (37) 
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dl]LW 


+  tR  [ALWRM2tw  +  (ALW  +  A  ALW)  RMHw) 

:|.  (2u)lW  ~  7SlwUlW  +  ^lw&lw)  /gg' 

(^lw  ~  2 8lw)2 

kDALWUJLw  {^LW  -  fow)  r  T  ,  .  X  r  .  ,  OA  B  'll 

- - -  [Ji(Aw)  {  -w  sin  r]L\v  +  2A^l  cos?]LW)J 

.  kDALWU!LW  {(jJlw  +  ^w)  (  T  r  a  \  ■  >  tt  r  a  \ 

H - - -  w  {  —  J\  (  Alw)  sm  rjLW  +  Hi{ALW)  cos  rjLW} 


+  2Azf  {J\{Alw)  cos  1JLW  +  Hi(Alw)  sin  r)LW} 
kn^Lw  {^lw  +  &lw)  9AAlw  f  WP 


driLW 


Vcp  (Alw  +  A  Alw) 


+  w  {J\  (. Alw  +  A Alw)  cos  t]lw  ~  Hi  (. Alw  +  A ALW)  sin  r/LW} 

kn^LW  {uLW  +  &Lw)  (A  .  A  A  \\  Wpd^A-LW 


(Alw  +  AAlw)  ycpJ 


drjLw 


f  dJ1(ALW  +  AAlw) 

w  | - drj~, - COS  ^ LW  ~  "'1  ^LW  +  A±Alw)  sm  rjLW 


9Hi(Alw  +  A  Alw) 


'  I  J_/W 

orjLW 

^kr,ujLw  (w lw  +  &lw)  9AAlw 

8  drjLW 

+  Hi  (Alw  +  A Alw)  cos  r]LW  1 


sin  r]LW  —  Hi  (. Alw  +  AAlw)  cos  tjlw 


2A zf  <  Ji  {Alw  +  A Alw)  svar)Lw 


kp^Lw  (^lw  +  &lw)  ,A  ,  a  A  \  To  A  b  f  9Ji(Alw  +  AAlw) 


(■ Alw  +  AAlw )  2AzL 


driLW 


sm  r]LW 


,  T  ,  A  ,  A  /»  A  I  QH^Alw  +  AAlw) 

+Ji  ( Alw  +  AAlw)  cos  t]lw  H - b - cos  tjlw 

ovlw 

—  Hi  (Alw  +  A ALw)  sin  rjLw  \ 
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dM, 


-B 

xLW 


dr)Lw-i 


k'D^LW  LW  +  ° LW )  9AAl\y 

8  dr)Lw-  i 


V^pP  (Alw  +  AAlw ) 


+  {</i  (-<4lw  +  AAlw )  cos  —  #1  {ARw  +  AAlw )  sin  txlvf} 


ko^LW  (a>lw  +  Cnv) 


(^Lty  +  AAlw ) 


cfcp 


dAA 


LW 


drjLw-i 


dJi(ALW  +  AAlw)  dHi(ALW  +  AAlw)  . 

w  < - - - cos  7?w - ~ - sin  ?7LW 

ot]lw- i 


drjLw-i 
ko^Lw  (^lw  +  Clw)  dAAi\y 
8  dr)LW_  i 

+  Hi  (Alw  +  AAlw)  cos 

.  kr)U>Lw  {alw  +  Cuv) 


2Az£  s  -A  (Alw  +  AAlw )  sin  77^ 


(^-w  +  AAlw) 


2Az 


b  i  9Ji(ALw  +  AAlw) 


dr/Lw-i 


dHi(ALW  +  AAlw ) 

H - a - cos 

oijlw- i 

E.  Pitching  Moment  Control  Derivatives 

The  pitching  moment  control  derivatives  are 

-B 


dMVRW  _  ArwUrw_PMi  +  (  /3 M2w  +  Vcp  KL  pm_ 


y%ph 


ddfiw  2 

-,B 


&ZRkL  I  a  2 


2 


3rw 


Arw  +  +  AARH/ 


a; 


w 


—  ^rw) 


where 


sm 

(40) 


(41) 


= 


(cos ax^ppkL  +  (sin ax}})®  +  Ax®}  /cd)  Ji(Arw)  cos r)RW 
~y^PFkLJi{ARW)  sinr]RW  — 


A->B 

^kj.A 


L^RW 


PM2rw  =  (cos ax™pPkL  +  (sin ax}})®  +  Ax®}  kD) 


A2RW  ^  -1)  (Arw)  cos  rjRW 


+  Hi(ARw)  sin  rjRw  \  +  ARw  ( ARw  +  AARw)  s  J\{ARw  +  AARw)  cos  qRw 


(42) 


—  H\(ARw  +  A ARw)  sin  qRw 


(43) 


P M3rw  —  A2mv  {Ji(Arw)  sin r]Rw  —  Hi(Arw)  cos  ijRw}  +  ARW  ( ARW  +  A ARW) 

* i -li  {Arw  +  A Arw)  sin  i]RW  +  Hi(Arw  +  A ARW)  cos  i]Rw \  (44) 
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dM, 


VRW 


du 


-Arw  (2 WRW  -  $Rw)_pM^  +  I  P M2RW  +  Vcp  thpM: 


RW 


4A 


RW 


4A 


RW 


3  RW 


Arw  +  Aw  +  AAwf))  f2^11'  7liR"'a,%  +  ^m3w  )(45) 


dA'L 


VRW 


-Arw^RW  pRW  —  Srw) 


dl]RW 


P rw  ~  2 5rw)2 

(cos  ax^pkL  +  {sin  ai^p  +  Ax kD) 


*  Jl(Arw)  sin  t)rw  —  ycp  kLJi(ARW)  cos tjrw 
+  — —  ^~~7~  +  ^  (cos ax^pPkL  +  (sin cxx^p  +  AxR}  A:d) 


4A 


.RW 


A 


{  —  Ji(Arw)  Anr]RW  +  Hi(Arw )  cos tjrw} 
8AArw 


rw- 


dijRw 


{Ji  ( Arw  +  AArW)  cos  rjRw  —  Hi  ( ArW  +  A ARW)  sin  t]rw } 


+  Arw  ( Arw  +  A  Arw)  |  —  J\(Arw  +  A  Arw)  sin  t]rw 
—  Hi{Arw  +  A  Arw)  cos  t/rw 


dJi  ( Arw  +  A  Arw) 


di]Rw 


cos  rjnw  - 


DH\  ( Arw  +  A  Arw  ) 

di]Rw 


sm  t)rw 


Vcp  k L^rw  pRW  +  &Rw) 

4  Arw 
8A  Arw 


A2rw  {Ji{Arw)  cost/rw  +  H\(Arw)  sin t)rw } 


4-Arw—z - {Ji  ( Arw  +  A  Arw)  sin  t]rw  +  Hi  ( Arw  +  A  Arw)  cos  i]rw} 

di]Rw 


+  Arw  ( Arw  +  AARW)  <^Ji(Arw  +  AARW)  cos  i]RW 
—  Hi(Arw  +  A  Arw)  sinr/ RW 

dJ\  (Arw  +  A  Arw)  .  dH\  (Arw  +  A  Arw) 

- sm  t]rW  H - - - cos  r)RW 


9i]rw 

AzRkLX>Rw  Prw  +  ctrvk) 


(Arw  +  A  Arw) 


dr]Rw 

dA  Arw 


drjRw 


(46) 
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8M, 


-B 

y  rw 


urw  {urw  +  &rw) 


dr/Rw- 1 


DM., 


-B 

y  lw 


A 


4A  RW 
8AARw 

RW- 


cos  ax^ppkL  +  {sino;a;^p  +  AxR  } 


{■h  ( Arw  +  AARW)  cos  r]RW  —  Hi  ( Arw  +  AARW)  sin  rjRW} 


dijRw-i 
+  ARw  ( ARw  +  A ARw) 

dJi  ( ARw  +  A ARw)  0H\  ( ARw  +  A ARw)  . 

*  < - - - cos  r]RW - - - sin  r)RW 


+ 


dr)Rw- 1 
yZP^L^RW  {^RW  +  Vrw) 


drjRw-i 


4  A 


RW 


+ARw t; - ■  {-h  (Aflvn  +  A ARw)  sin  r)Rw  +  Hi  {ARw  +  AAr^v)  cos  i]rw} 

orjRw- i 

+  ArvF  (ArvK  +  AArw 

<9Jl  (Arw  +  AArw) 


dr)Rw- 1 


dH\  (Arw  +  A Arw) 
sin  l/rw  H - - - cos  r]RW 


A  ZR  kLUJRw  {x>RW  +  CTRiy) 


(Arw  +  AArw 


drinw-i 

8  A  ARw 


85 


ALWULWpM  +  f 


w 


4A 


w 


4A 


lw 


drjnw-i 

3lw 


H - ^  (  A|w  +  (Arw  +  A Alw)2 


u 


LW 


(^lw  ~  25lwY 


where 


(47) 


(48) 


= 


(cos  ax^pkL  +  {sinax^p  +  Axf  }  kD )  J\(ALW)  cos r]LW 

Azf 

~y^pF YlJi{Alw)  sin  ijlw  - 


^L-kr,A 


LSiLW 


PM2r 


=  (cos  ax^PkL  +  {smax^p  +  Axf  }  /cr) 


Afw  l  Ji(Alw)  costjlw 


+  Hi(Alw)  siwrjLw  >  +  Arw  (Arw  +  A ALW)  <  Ji(Alw  +  AALW)  cos rjLW 


(49) 


—  H\{Alw  +  AAlw)  sin  ??rw 


(50) 


P Mzlw  —  A2lw  {  J\(Alw)  sin  r/LW  —  Hi(Alw)  cos  l/lw}  +  ALW  ( ALW  +  A ALW) 

Ji(Alw  +  A  Alw)  sin  t/lw  +  Hi(Alw  +  A  Alw)  cos  t/rw  1  (51) 
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dh'L 


-B 

VLW 


—ALw  {2^lw  —  5lw) 


du> 


P  i.w  + 


PM2lw 


LW 


-B 


4  A 


LW 


4  A 


LW 


3lw 


A-ZftkL  (  A2  ,  A  ,  A  A  f^ULW  7&lwu\w  +  ^LW^LW 

—  (Alw  +  (ALW  +  AAlw)  )]  ^  —  -  —2 

cos  ax^ppkL  +  {sin  ax^pp  +  Axp }  kD) 


(52) 


8MP lw  _  —Alwujlw  (dlw  —  8Lw) 
drjLw  2 

*  J\{Alw)  sin rjiw  —  y^pPkLJi{ALw)  cos  tilw 

.  ^lw  {^lw  +  &lw) 


4  A 


LW 


(cos ax^p kL  +  {sin ax^p  +  Ax^}  ko) 


+  A 


A2lw  {  —  Ji(Alw)  sin  rjLW  +  Hi(Alw)  cos  rjLW} 
8AAlw 


lw~ 


8t]lw 


{Ji  (Alw  +  A Alw)  cos  qLw  —  H i  (Alw  +  A Alw)  sin  i]lw} 


+  Alw  ( Alw  +  A Alw)  |  —  Ji(Alw  +  A Alw)  sin rjLW 

—  Hi(Alw  +  A  Alw)  cos  i]lw 


dJi  ( Alw  +  A Alw) 


drjLw 


cos  7]LW  - 


8H\_  ( Alw  +  A Alw) 

dl]LW 


sin  r)LW 


VcpPkLU)LW  (ulw  +  &lw) 


+Alw 


4  ALw 
8AAlw 


A2lw  {Ji(Alw)  cos  rjLW  +  Hi(Alw)  sin  t]lw} 


drjLw 


{J\  ( Alw  +  A  Alw)  sin  rjLW  +  Hi  (. ALW  +  A  ALW)  cos  r]LW} 


+  Alw  ( Alw  +  A Alw)  l  J\  {Alw  +  A Alw)  cos  qLw 

—  Hi{Alw  +  A Alw)  sin r)LW 

8J\  ( Alw  +  A  Alw)  .  8H\  {Alw  +  A  Alw) 

- sin  rjLW  H - ~ - cos  rjLw 


8t]lw 

AzPkLLJLW  (ULW  +  &Lw) 


( Alw  +  AAlw) 


8l]lw 
8  AAlw 


8t)lw 


(53) 
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dM, 


-B 

y  lw 


UlW  \WlW  +  a LW ) 


drjLw-i 


A 


4  Alw 

BAAlw 

LW' 


cos  ax^ppkL  +  {sin  ax^p  +  AxR]  ko) 


{■h  ( Alw  +  AAlw )  cos  t)lw  ~  Hi  ( ALW  +  A ALW)  sin  rjLW} 


drjLW-i 

+  Alw  ( Alw  +  A Alw) 

(  dJi  ( Alw  +  A  Alw)  dHi  ( Alw  +  AAlw )  . 

*  <  - « - cos  r]LW - — - sm  r)Lw 

L  orjLw- 1  ot]lw-  1 

,  y^PkL^LW  (^LW  +  &Lw) 


+Alw 


4:ALW 

dA  Alw 


driLw-i 

+  ALw  (Alw  +  AALw) 

dJ\  (Alw  +  A  Alw 


{  J\  (Alw  +  AAlw)  sin  t)lw  +  Hi  (ALw  +  A ALw)  cos  t]lw} 


drjLw-i 

AZftkLUJLW  (ulw  +  o~lw) 


dHi  (Alw  +  A  Alw) 
sm  rjLw  H - 5 - cos  rjLw 


(Alw  +  AAlw 


drjLw- 1 
dA  Alw 


drjLw- 1 


(54) 


F.  Yawing  Moment  Control  Derivatives 


dM. 


-B 

zRW 


dd 


RW 


- ^ — -^-Y M\  Rw  +  -  [koARWY  H2rw  —  kLARWY Mzrw 

+  ko  (Arw  +  A Arw)  Y MArw  —  kL  (ARw  +  A ARW)  Y M5rw\ 


lo 


RW 


dM. 


zRW 


Xurw  ~  2 5rw)2 
Arw  (2 (jJrw  ~  &rw 


(55) 


dco 


RW 


-Y Mirw  +  -  (kDARWY M2rw  —  kLARWY M3rw 


+  kD  (Arw  +  A Arw)  YMArw  —  ki  (Arw  +  AArw)  Y M5rw) 

2cU RW  ~  7§R\y&Rw  A  4c iJrw^rw 

(urw  —  2  5rwY 


(56) 


Approved  for  public  release;  distribution  unlimited. 
17  of  36 


where 


ym1rw 

ym2rw 

ym3rw 

ymArw 

ym5rw 


Ji(Arw)  sin tjrw  \  —kD  (sinaa;^p  +  Axp)  —  cos  ax^ppkL 


vWPk 
Ucp  A 


cos  r)RW 
'  sin  r)RW 


_  wkiAmv 
4 

(sinaxjp  +  AxR)  {  Ji{Arw )  sin r]RW  -  #i(Aw)  cosrjRW} 
-cos axJp{Ji(Aw)sin77W  -  Ih(ARW)  cosr]RW} 

+y^P  {Ji(Arw)  cosijrw  +  Hi(Arw)  sin  r/RW/}  H - 

(sin  ax™p  +  Ax|)  |  Ji(Arw  +  AARW)  sin 


(57) 

(58) 

(59) 


+  Hi(Arw  +  AArw)  cos rjRW  j 

—  cos  ax^pp  {  J\{Arw  +  A Arw)  sin  i]RW  +  H\  (Amv  +  AARW )  cos  r/Rw} 

+?/cpP  {•MArvp  +  ^Arw)  cosrjRW  —  Hi{Arw  +  A ARW)  sin r/Rw} 
w  {Arw  +  A Arw} 

4 


(60) 


(61) 
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-B 


9MzRW  _  Arw^rw  (urw  —  Srw) 
drinw  2 


Ji(Arw)  cosri{—kD  (sinax^p  +  Ax^) 


-  cos aix™pkL}  +  Ji(ARW)y™pkLsinr)RW 
.  koARw<jJRw  ( ujrw  +  o~Rw) 


sm  ax 


WP 


4  v - cp 

*  {J\{Arw)  cos r]Rw  +  H\(ARw )  siniy^^) 
kL^RW^RW  rw  +  prw) 


+  Ax|) 


cos ax^p  {■J\{Arw)  cos rjRW  +  Hi(Arw)  sin r/RW} 


+  Vcp  {  —  Ji(Arw)  sin ijrw  +  Hi(Arw)  cost/rw} 
kD^RW  (urw  +  <trw)  (sin  ax^p  +  Ax|) 


8AA 


RW 


dr/nw 


{Ji{Arw  +  A Arw)  sin  r//?vF  +  Hi(Arw  +  A ARW)  cos rjRw} 


+ 


{Arw  +  AARW )  <j  - — — — - - - -  sin  r/RW  +  Ji(Arw  +  AAKM/)  cos  rjRw 

OVRW 


dHi(ARw  +  AArw) 
dr)RW 

kL^Rwi^RW  +  ctrw) 


cos  rjRw  —  H\  {Arw  +  AArw  )  sin  r/pvv 

9AARw 
di]Rw 


—  cos  axcp  (  J\  {Arw  +  AARW)  sin  r/RW 


H\{ARw  +  A ARw)  cos  i]rw) 


+  2/Cp  Ji(Arw  +  AARW/')  cost/rw  —  H\  (ARW  +  AApiy)  sin  r/pyv 

+ 

+ 


w{Arw  +  A  Arw) 


,A  i  a  4  \f  h/p  ( dJi(ARW  +  AArw)  . 

[Arw  +  A  A  w )  I  -  cos  axcp  I  - — — - sm  r/RW 

,  .  .  .  .  d HA Arw  +  AApjy) 

Ji{Arw  +  AApvi/)  cos ?7pvk  H - ^ - - —  cos  t/rw 


dr/mv 


+  y , 


H\{Arw  +  AArW)  AnrjRw 
wp  i  dJi(ARW  +  A  Arw 


cp 


dyRw 


cos  rjRw  —  -h  {Arw  +  AArw)  sin  rjRw 


d Hi( Arw  +  AArw)  .  u  t  a  ,  A  a  ^  \  wdAARW 

- sm  r/pvi'  —  Hi\Arw  +  AARW)  cos  rjRW  + 


dr/Rw 


4  cfyw 


(62) 
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9MzRW  kDURw  {wrw  +  ctrw)  (sin  ax^pP  +  Ax  ft) 
drjRw-i  4 


8A  A 


RW 


{J\(Arw  +  AArw )  sin  rjRw  +  H\  (ARW  +  A ARW)  cosriRW} 
dJ\{ARw  +  A  ARw) 


drjRw- 1 


drjRW- 1 

+  (Arw  +  AArw ) 

dHi(ARw  +  AArw) 

H - o - cos  rinw 

oVrw-i 

_  ^L^Rwi^RW  +  Orw)  <9AArw 
4  _dr)Rw-i 

+  Hi(Arw  +  A Arw)  COS  ??rw) 


sin  r/Rw 


—  cosaxcp  (  Ji(Arw  +  AArw)  sin  rjnw 


+  2/cpP  Ji(Arw  +  AArw)  cosrjRW  —  H\  {ARW  +  A ARW)  sin  rjRw 


?  CP 


+ 


w(Arw  +  A  Arw) 


+  (Arw  +  AARw)  |  —  cos  axcp 
dHi(ARW  +  A  Arw) 


wp  ( dJi(ARW  +  A ARw) 


drjRw-  i 


sin  riRW 


9rjnw-\ 


COS  TjRW 


y , 


wp  /  dJi(ARW  +  A  Arw) 


cp  \  dyRW-i 
dHi(ARW  +  A  Arw) 
drjpw- 1 


COS  f]RW 

w  dAARw 

4  drjnw-i 


sm  rjRW  + 


(63) 


For  the  left  wing,  the  yawing  moment  control  derivatives  become 

-B 


8M. 


zLW 


Alw^lw 


YM1 


1 

LVF  ^ 


IzdAlwY  M2lw  —  kLALWy  M3lw 


9Slw  2 

+  kn  (Alw  +  AALW)  Y MAlw  —  kL  ( ALW  +  A ALW)  Y M5lw 


u 


LW 


(^lw  —  25 lwY 


(64) 


-B 


9MzLW  _  —Alw  (2cn lw  —  ^LW^YM  —  - 
Oujlw  2  LW  4 


kpAm/Y  M2lw  —  UlAlwY  M3lw 


+  kR  (Alw  +  AAlw)  y  M±lw  —  kL  ( ALw  +  AALw)  i'  M5lw 

2 ^lw  —  75lw^lW  4“  4 ujlw^lw 
(^lw  ~  2 5rw) 


(65) 
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where 


ym2lw 

ym3lw 

YMiLW 

YMbLW 


J\{ALW)  sin  t]LW  <  —  kD  (sincur^p  +  Axf)  —  cos ax^ppkL 


vWPh 
y cp  '  ■ 


cos  r]LW 
'  sin  rjLW 


wUlAlw 

4  ^ 

(sin  ax^p  +  Axf)  {  Ji(Alw )  sin  rjLW  -  H^Alw)  cos  rjLW} 
-  cos  ax^vP  {Ji{Alw)  sin  rjLW  -  #i(Aw)  cos  rjLW} 

+y^>P  {^i (Alw)  cos r]LW  +  Hi(Alw )  sin77W}  H - -AAL 

(sin  ax^p  +  Axf )  j  Ji{Alw  +  AALW )  sin  rjLW 


(66) 

(67) 

(68) 


H\{Alw  +  AAlw )  cos tilw 


(69) 


—  cos  ax^pP  {Ji(Alw  +  AAlw)  sim]LW  +  Hi(Alw  +  AALW )  cosr/ivy} 

Ay^pP  {Ji(Alw  +  AALW)  cosjjlw  —  Hi(Alw  +  AALW)  sinry^} 
w  (Any  +  AAj^y} 

4 


(70) 
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dMzLW  Alwujlw  ( uLw  -  5lw )  T  ,  .  ,  (if-  wp  ,  a  b\ 

= - - -  J\  [A lw  )  cos  77  (  -  kD  (sin  axcp  +  AxL ) 


-  cos  ax^pkL}  +  Ji(ALW)yWp kLsmr)LW 


kr> Alwujlw  {ulw  +  °lw)  /  .  wp  ,  A  b\ 
- - - (sm  axcp  +Axl) 


*  (  J\{Alw)  cos t]lw  +  Hi(Alw)  si nrjLw 


-  cos  axZP  {Ji(Alw)  cos r)LW  +  Hi(Alw)  sin rjLW} 


JclAlW^LW  (ulW  +  &  LW ) 


+  V^vF  {— Ji(Alw)  sin  rjLW  +  H\  (ALW)  cos  t)lw}  J 

kDULW  (uLw  +  &lw)  (sin  ax™P  +  Axf ) 

4 

*  — q — ~  {-MAlw  +  AAlw)  simjLw  +  Hi(Alw  +  AALW)  cosr]Lw} 
orjLw 

,  ^  f  dJi(ALW  +  AAlw)  .  ,  w ,,  ,  A  /,  ^ 

+  (  Alw  +  AAlw)  < - - sm rjLW  +  Ji(Alw  +  AALW )  cos r//-jW 


8Hi(Alw  +  AAlw ) 


ovlw 

kL^Lw{uLW  +  0~w)  OAAlw 
4  [ 

+  H\{ALw  +  AAlw)  cos  r)Lw) 


cos  t]lw  ~  H\{ALw  +  AALw)  sin  tjlw 


-  cos  a£p  (Ji(Alw  +  AAlw)  sin r)Lw 


+  ?/«^P  |  -h  (ALw  +  AAlw)  cos  t]lw  —  Hi{ALw  +  AALw)  sin  ijLw 
w(ALw  +  AAlw)  \ 


+  ( ALw  +  AAlw)  s  -  cos  ax: 


wp  ( 8Ji(ALw  +  AAlw) 


dr]Lw 


sm  TjLw 


T  i  \  ,  a  /i  \  ,  8Hi(ALw  +  AAlw) 

+  Ji\Alw  +  AAlw)  cos  tjlw  H - 77 - cos  rjLw 

dr/LW 

—  Hi(Alw  +  AAlw)  sin  i]lw  ^ 

WP  f  dJi(ALw  +  AAlw)  t,a  .  A  ,  \  • 

+  yCp  I - ~drj— - cos 7 ~ilw  ~  Ji\ALw  +  AALw)  sm rjLw 

dHi(ALw  +  AAlw)  .  u  t  a  ,  A  a  \  \  iv  dAALw 

—  - o - sin —  H\\  Alw  +  AALw)  cos  qLw  +  -7—^ - 

J  4  &qLW 
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dM, 


-B 

zLW 


kDuLW  ( uLW  +  aLW )  (sin  ax^p  +  Axf ) 


drjLw- 1 


<9AA 


{A(Ai,w  +  AA^vy)  sin tjlw  +  Hi(Alw  +  AApjy)  cos  r/ 7-4,5/} 
9Ji(Alw  +  AApyy) 


dr)Lw- 1 


dr/Lw- 1 

+  (Alvk  +  AAlw) 

OHi(Alw  +  AApyy) 

H - o - cos 

OljLW- 1 

kL^Lwi^LW  +  OlhO  ^AA^yy 
4 

+  Hi(Alw  +  AAlw)  cos  ??lw) 


sm  ?/Lw 


-  cos  cnc^p  (Ji(Aw  +  AAlw)  smrjLW 


+  ?/™P  Ji(Alw  +  AAlw)  cost/lw  —  Hi(Alw  +  AAiM/)  sinr/pyy 


(72) 


w(Alw  +  AAlw) 


+  (Aw  +  AAlw)  j  -  cos  ai^ 
9Hi(Alw  +  AAlw ) 


vup  / dJ\(ALw  +  AApW) 


drjLw- 1 


sm  r/Mv 


dljLW-l 


cos 


+  V, 


WP  f  dJi(ALw  +  AApvi/) 


cp  V  dr/Lw- 1 

9Hi(ALw  +  A  Apy^) 
dVLw-i 


cos 

ic  DAAlw 
4  dr/Lw-i 


smr/Lw  + 


III.  Aerodynamic  Control  Derivatives  about  Hover 

The  control  of  this  vehicle  in  the  vicinity  of  hover  is  of  considerable  interest.  Therefore, 
the  control  derivatives  are  evaluated  at  the  hover  condition  where  =  culiv  =  cu0,  = 
=  0,  ?7pyi/  =  r/pvi/  =  0,  and  AArw  =  A^Alw  =  0.  Additionally,  it  is  assumed  that 
the  nominal  center-of-gravity  of  the  vehicle  and  wing  root  hinges  are  aligned  such  that 
Azbr  =  A zB  =  0.  When  the  nominal  center-of-gravity  and  wing  root  hinges  are  not  aligned,  a 
non-zero  cycle-averaged  pitching  moment  is  produced  which  rotates  the  lift  vector,  translates 
the  vehicle,  and  does  not  allow  hover.  Lastly,  since  A  A  =  77  —  r/_i,  =  1  and  |AT  —  —  1. 
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A.  X-Body  Axis  Force  Control  Derivatives  About  Hover 


dFxRw 

85rw 

9FxLW 

95lw 

dFxRW 

dujRw 

9FxLW 

Oujlw 

®FxRW 

dr]Rw 

9FxLW 

drjLW 

9FxRW 
drjRw- 1 

dFxLW 
drjLW- 1 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


kL^-LW^o 


k-^A, 


k-^A, 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 
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B.  Y-Body  Axis  Force  Control  Derivatives  About  Hover 


dF 

U1  v  rw 

9Srw 

ggL 

ddLw 

dFymv 

dujRw 

®FyLw 

Oivlw 

dFymv 

dijRw 

ggL 

dl]LW 

dFyRw 

drjRw-i 

dFyhw 

d^LW-l 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


o 

o 

0 

0 


kDul 


4 

kDu2Q 

4 

kD^o 


4 

knujl 


Hi(Arw) 
Hi(Alw ) 
Hi(Arw) 


Hi(Alw) 


(81) 

(82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 
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C.  Z-Body  Axis  Force  Control  Derivatives  About  Hover 


dFzRw 

35rw 

dFzLw 

B5lw 

3Fzrw 

dujRw 

®Fzlw 

duJLW 

®Fzrw 

dr/Rw 

9FzLW 

drjLW 

®FzRW 

driRW- 1 

&FzLW 

d^LW-l 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


— kp  ArW  J i  (  Arw  )  uj0 
~kr)ALwJi{ALw)^o 
0 

0 

—kpJ\{ARw)^ 

4 

~ kpJi(ALW)ujl 
4 

kpJi(ARW)ujl 

4 

kpJi(ALW)ujl 

4 


(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 
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D.  Rolling  Moment  Control  Derivatives  About  Hover 


dM 


-B 

XRW 


95rw 

b 


dM, 


xLW 


95lw 

B 


dM, 


XRW 


dujRw 

-B 


dM ■ 


xLW 


Qujlw 

B 


dM, 


XRW 


dr]Rw 
9MX  LW 


dr)Lw 

B 


dM 


xRW 


dr)Rw- 1 
dMx  LW 


drjLw-i 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


hover 


__d_^u — o  | yWpARW  +  wJi(Arw)] 
D  7p - °  \VqpP A-LW  +  wJ\{ALW )] 

0 

0 

— [2 yWpARW  +  wJi(Arw )] 

~  [2  yWpALw  +  wJi(Alw)] 

~  [2^PAijvr  +  wJi{Arw)\ 

— ^r2-  [2 yWpALW  +  wJi(Aivv)] 


(97) 

(98) 

(99) 
(100) 
(101) 
(102) 

(103) 

(104) 
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E.  Pitching  Moment  Control  Derivatives  About  Hover 


=  ATiWJx(ARW)uj0  (cos ax^pkL  +  (sin ax^p  +  Axf  }  kD)  (105) 
=  AlwJi(Alw)u)0  (cos ax^pkL  +  (sinax^p  +  Axf  }  kD)  (106) 


_ lv±VRW 

du>Rw 


(107) 


-,lv±yLW 

diOLW 


(108) 


drjRw 


ARwulyWpkLJi{ARW)  + 


Jlvly  lw 

dl]LW 


+  y^pP  ^lHi{ArW) 


ALw^yWp  kLJi(ALW)  + 


+  y^>P ^lH\{Alw) 


drjRw-i 


y^ppkLHi(ARW) 


driLw-i 


+  y^>P kLH\{ALW) 


(109) 


*  ^cosax^pkL  +  (sin cxx^pp  +  Ar^}  kD  ^  Ji(Arw)  (110) 


(111) 


(112) 


*  ^cos  ax^p  kL  +  {sin  ax^p  +  Axf}  kD^j  Ji(Alw)  (113) 


(114) 


(115) 


*  (^cosax^pkL  +  {sinox^p  +  A;rp}  kD  ^  Ji(Arw)  (116) 


(117) 


(118) 


*  ^cos  ax™pPkL  +  (sin  ax^pp  +  Axf }  kD  ^  J\(ALW)  (119) 


(120) 
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F.  Yawing  Moment  Control  Derivatives  About  Hover 


9MzRW 

OSrw 

hover 

9MzLW 

35lw 

hover 

9MzRW 
Blorw  , 

hover 

9MzLw 

8ujlw 

hover 

9MzRW 

hover 

9MzLW 

dl]LW  , 

hover 

9MzRW 
drjRw-i  , 

hover 

9MzLW 

dl]LW-l 

hover 


o 

o 

—2 kiu0  ARwyWpJl(Anw)-* - p — 

2k.LU0  ALwy'^P Ji(Anv)  H - p — 

ppffi (Amv)  (smax"F  +  Arrg) 

~(~  cos  ax?/ H^A^)  +  vZPMA«w)  + 

-  H|M/,ni  (sin  axj^F  +  Aif) 
-*%£(-«*  axY/H,  (Alw)  +  y%F.h  (ALW)  + 

-  H i !  A Itw  i  (sin  nij'’  +  Ai(() 

(- cos ax^p Hi {Ayw)  +  vI-UAkw)  + 
H !  i  4  in  i  (sin  QX^P  +  Ai£) 

-  ^  COS  (diB.)  +  4  P  J,  (21„„)  + 


(121) 

(122) 

(123) 

(124) 

(125) 

(126) 

(127) 

(128) 

(129) 

(130) 

(131) 

(132) 
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G.  Control  Effectiveness  Matrix 

Using  the  general  form  of  Equation  2,  writing  this  expression  for  all  three  forces  and  moments, 
and  separating  out  the  control  and  non-control  parameters  yields 


$rw 

&LW 

A  F*des 
am: is 

—  Bai 

A  (jJrw 

Ac olw 

a  m 

Vrw 

AMzBdes_ 

.  Vlw  . 

~VRW-1 

~Vlw-i 


(133) 


where  Bai  is  the  control  effectiveness  matrix,  that  contains  the  aerodynamic  control  deriva¬ 
tives  evaluated  at  hover,  and  the  vector  on  the  left  side  of  Equation  133  is  a  desired  set 
of  forces  and  moments,  which  are  generated  by  a  cycle-averaged  control  law  that  will  be 
discussed  shortly,  to  achieve  a  desired  maneuver.  Since  Vrw-IiVlw-\  are  not  control  vari¬ 
ables,  B A2  is  not  considered  a  control  effectiveness  matrix.  Instead,  the  terms  containing 
Vrw-u  Vlw-i  simply  produce  a  set  of  forces  and  moments  that  are  only  persistent  over  a 
single  wingbeat  cycle.  Hence,  Equation  133  can  be  manipulated  to  give 


'A  Hi,' 

$RW 

A  F^ls 

5lw 

A  FzBdes 
AKls 

—  B  A2 

~Vrw-i 

—IjLW-l 

—  Bai 

Ac orw 

Ac olw 

A  i\  r  B 

Vrw 

AMll,. 

Vlw 

(134) 
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The  matrices  Bai  and  Ba2  are  explicitly  expressed  as 
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OSrw 
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0 
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From  the  structure  of  Bai,  it  is  found  that  this 


matrix  has  rank  six. 


Hence,  split-cycle 


constant  period  frequency  modulation  with  wing  bias  has  provided  a  technique  to  control 
6  degrees-of- freedom  with  only  two  physical  actuators.  Fundamental  frequency  (ujrw^lw) 
can  be  used  to  control  the  x-body  axis  force  and  yawing  moment.  Likewise,  the  split-cycle 
parameters  (8rw-,Slw)  can  be  used  to  control  rolling  moment  and  either  z-body  axis  force 
or  pitching  moment.  The  wing  bias  parameters’  primary  function  is  to  replace  the  bob- 


weight3,4  as  a  pitching  moment  control  parameter.  The  control  allocation  objective  is  to  find 

1  T 


the  control  vector, 


Srw  $lw  Aljrw  Aulw  tIrw  Vlw 


,  such  that  Equation  134  is 


satisfied.  The  desired  forces  and  moments  in  Equation  134  are  generated  by  a  feedback 


structure  designed  to  track  a  trajectory.  The  control  allocation  scheme  is  a  pseudo- inverse 
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and  the  control  vector  is  computed  as 


'a  Fils' 

\ 
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(137) 


where  Bai+  is  the  pseudo-inverse  of  Bai  The  control  values,  in  Equation  137,  are  held 
constant  over  a  wingbeat  cycle  and  are  applied  to  the  oscillator  generating  the  wingbeat 
forcing  function. 


IV.  Results 


The  modeling  and  control  analysis  is  now  utilized  in  a  six  degree-of- freedom  (DOF) 
model  of  the  flapping  wing  MAV.  The  control  design  approach  is  the  same  as  that  used  by 
Oppenheimer,  et.  al.4  Figure  1  shows  the  structure  of  the  controller.  Position  or  attitude 
errors  generate  desired  accelerations  (or  forces  and  moments).  The  forcing  functions  are 


Figure  1.  Five  degree-of- freedom  split-cycle  controller  with  wing  bias. 

inputs  to  an  aerodynamic  model  of  the  vehicle  and  6  DOF  equations  of  motion  are  integrated 
to  obtain  the  state  of  the  vehicle.  The  6  DOF  model  is  driven  by  the  instantaneous  forces  and 
moments,1  from  the  blade-element  model,  and  not  the  cycle- averaged  forces  and  moments. 
The  commands,  for  the  presented  simulation  run,  are  shown  in  Figure  2.  The  vehicle 
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Figure  2.  Commanded  MAV  Trajectory. 

starts  by  executing  a  roll  to  align  its  heading  with  waypoint  1  and  then  pitches  over  to 
perform  a  translation  in  the  direction  of  the  y-axis  of  the  inertial  frame.  Following  this, 
the  vehicle  again  aligns  its  heading  with  waypoint  2  and  proceeds  to  translate  in  both  the 
x-axis  and  z-axis  directions  of  the  inertial  frame.  Two  constant  altitude  maneuvers  are  made 
followed  by  a  change  in  altitude  to  get  to  the  final  waypoint  (5).  Figure  3  shows  the  roll, 
pitch,  and  yaw  angles  of  the  MAV  during  the  flight,  while  Figure  4  shows  the  Euclidean 
norm  of  the  position  errors.  The  large  roll  angles  are  used  to  align  the  heading  of  the  vehicle 
with  the  next  waypoint,  while,  for  forward  translation,  a  non-zero  pitch  attitude  angle  is 
utilized.  The  position  errors  are  small  with  a  maximum  of  about  0.12  m. 

Figures  5,  6,  and  7  show  the  values  of  the  control  variables  required  to  produce  this 
maneuver.  The  trim  wingbeat  frequency  is  about  120  Hz  with  large  excursions  during  the 
portions  of  the  flight  where  altitude  changes  are  made.  The  split  cycle  parameter  is  used  to 
produce  roll  and  is  therefore  active  during  all  portions  of  the  flight  except  movement  to  the 
last  waypoint  when  there  is  no  heading  error.  Similarly,  the  bias  is  used  to  pitch  the  vehicle 
to  assist  forward  translation.  Hence,  the  bias  is  active  during  all  portions  of  the  flight.  One 
interesting  point  here  is  that  a  bias  value  of  less  than  about  0.6°  is  sufficient  to  regulate  the 
attitude. 
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Euler  angles  vs.  time 


Figure  3.  Vehicle  Attitude. 


Tracking  error  vs.  time 


Figure  4.  Euclidean  Norm  of  Position  Error. 


V.  Conclusions 

In  this  work,  a  control  strategy  for  a  flapping  wing  micro  air  vehicle  was  developed.  The 
control  strategy  made  use  of  split-cycle  constant  period  frequency  modulation  with  wing  bias 
to  provide  six  DOF  control  over  the  vehicle.  Simulation  results  show  that  sufficient  control 
authority  exists  to  track  six  DOF  trajectories.  The  cycle-averaged  controller  provides  the 
ability  to  track  desired  trajectories  when  applied  to  a  blade-element  based  simulation  model 
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1000 


Fundamental  frequencies  vs.  time 


Figure  5.  Fundamental  Wingbeat  Frequency,  lo. 


Split  cycle  parameters  vs.  time 


Figure  6.  Split-Cycle  Parameter,  S. 


of  the  flapping  wing  MAV  where  the  instantaneous  aerodynamic  forces  and  moments  forced 
the  equations  of  motion. 


References 

1Oppenheimer,  M.  W.,  Doman,  D.  B.,  and  Sigthorsson,  D.  O.,  “Dynamics  and  Control  of  a  Biomimetic 
Vehicle  Using  Biased  Wingbeat  Forcing  Functions:  Part  I  -  Aerodynamic  Model,”  Submitted  to  2010  AIAA 

Approved  for  public  release;  distribution  unlimited. 

35  of  36 


Wing  biases  vs.  time 


Figure  7.  Wing  Bias,  77. 


Aerosciences  Conference,  2010. 

2Wood,  R.  J.,  “The  First  Takeoff  of  a  Biologically  Inspired  At-Scale  Robotic  Insect,”  IEEE  Transactions 
on  Robotics ,  Vol.  24,  No.  2,  2007,  pp.  341-347. 

3Doman,  D.  B.,  Oppenheimer,  M.  W.,  and  Sigtlrorsson,  D.  O.,  “Dynamics  and  Control  of  a  Minimally 
Actuated  Biomimetic  Vehicle:  Part  I  -  Aerodynamic  Model,”  AIAA  Guidance,  Navigation  and  Control 
Conference,  AIAA-2009-6160,  Aug.  2009. 

4Oppenheimer,  M.  W.,  Doman,  D.  B.,  and  Sigtlrorsson,  D.  O.,  “Dynamics  and  Control  of  a  Minimally 
Actuated  Biomimetic  Vehicle:  Part  II  -  Control,”  AIAA  Guidance,  Navigation  and  Control  Conference, 
AIAA-2009-6161,  Aug.  2009. 

5Gradshteyn,  I.S.  and  Ryzhik,  I.M.,  Table  of  Integrals,  Series  and  Products,  6  ed.,  Academic  Press,  2000. 


Approved  for  public  release;  distribution  unlimited. 

36  of  36 


